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Abstract The fractional calculus of variations is a subject under strong research. 
Different definitions for fractional derivatives and integrals are used, depending 
on the purpose under study. The fractional operators in this paper are defined 
in the sense of Jumarie. This allows us to work with functions which are non- 
differentiable. We present necessary and sufficient optimality conditions for frac- 
tional problems of the calculus of variations with a Lagrangian density depending 
on the free end-points. 



1 Introduction 



The Fractional Calculus (FC) is one of the most interdisciplinary fields of mathe- 
matics, with many applications in physics and engineering. The history of FC goes 
back more than three centuries, when in 1695 the derivative of order a = 1/2 was 
described by Leibniz. Since then, many different forms of fractional operators were 
introduced: the Grunwald-Letnikov, Riemann-Liouville, Riesz, and Caputo frac- 
tional derivatives EUEIESI), and the more recent notions of ifTTl [191 1241 |25l . FC 
is nowadays the realm of physicists and mathematicians, who investigate the useful- 
ness of such non-integer order derivatives and integrals in different areas of physics 
and mathematics (see, e.g., ifTUl [1811231 ). It is a successful tool for describing com- 
plex quantum field dynamical systems, dissipation, and long-range phenomena that 
cannot be well illustrated using ordinary differential and integral operators (see, e.g., 
[15, 18 24 36]). Applications of FC are found in classical and quantum mechanics, 
field theories, variational calculus, and optimal control (see, e.g., |[T4l[T7ll20ll ). 

The calculus of variations is an old branch of optimization theory that has many 
applications both in physics and geometry. Apart from a few examples known since 
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ancient times such as Queen Dido's problem (reported in The Aeneid by Virgil), 
the problem of finding optimal curves and surfaces has been posed first by physi- 
cists such as Newton, Huygens, and Galileo. Their contemporary mathematicians, 
starting with the Bernoulli brothers and Leibniz, followed by Euler and Lagrange, 
invented the calculus of variations of a functional in order to solve those prob- 
lems. Fractional Calculus of Variations (FCV) unifies the calculus of variations 
and the fractional calculus, by inserting fractional derivatives into the variational 
integrals. This occurs naturally in many problems of physics or mechanics, in or- 
der to provide more accurate models of physical phenomena. The FCV started 
in 1996 with the work of 11361 . Riewe formulated the problem of the calculus of 
variations with fractional derivatives and obtained the respective Euler-Lagrange 
equations, combining both conservative and nonconservative cases. Nowadays the 
FCV is a subject under strong research. Different definitions for fractional deriva- 
tives and integrals are used, depending on the purpose under study. Investigations 
cover problems depending on Riemann-Liouville fractional derivatives (see, e.g., 
@[[5][l6][32)), the Caputo fractional derivative (see, e.g., 0] EES] ED), the sym- 
metric fractional derivative (see, e.g., 11241 1. the Jumarie fractional derivative (see, 
e.g., HE] [E] III El] |22] |27)), and others EQZDCH). For applications of the frac- 
tional calculus of variations we refer the reader to |[T5l l20l l24l [351 - Although the 
literature of FCV is already vast, much remains to be done. 

In this paper we study problems of FCV which are defined in terms of the Ju- 
marie fractional derivatives and integrals. The Euler-Lagrange equations for such 
problems with and without constraints were recently shown in 0. Here we develop 
further the theory by proving necessary optimality conditions for more general prob- 
lems of FCV with a Lagrangian that may also depend on the unspecified end-points 
y(a), y{b). More precisely, the problem under our study: to extremize a functional 
which is defined in terms of the Jumarie fractional operators and having no con- 
straint on y(a) and/or y(b). The novelty is the dependence of the integrand L on the 
a priori unknown final values y(a), y(b). The new natural boundary conditions ©- 
(O have important implications in economics (see lfl2l and the references therein). 

The paper is organized as follows. Section E] presents the necessary definitions 
and concepts of Jumarie's fractional calculus. Our results are formulated, proved, 
and illustrated through examples in Section [3] Main results of the paper include 
necessary optimality conditions with the generalized natural boundary conditions 
(Theorem[T]i that become sufficient under appropriate convexity assumptions (The- 
oremE])- We finish with Section[4]of conclusions. 



2 Fractional Calculus 

For an introduction to the classical fractional calculus we refer the reader to 11231 
EU [34] [38) . In this section we briefly review the main notions and results from the 
recent fractional calculus proposed by Jumarie 119. 201 I2TI . 
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Definition 1. Let / : [a,b] — > M be a continuous function. The Jumarie fractional 
derivative of / is defined by 

f (a) V := TT-7T\ (\t-x)- a -\f{x)-f{a))dT, a<0, 
1 (-a) Jo 

where F(z) ~ t z ~ l e~' dt . For positive a, one will set 

/(«)(,) = (f^- l \t))' = Y^—jj t J\t-x)-<*(f(x)-f(a))dx, 
for < a < 1, and 

/ (a) (f):=(/ (a ~ n) W) W , n<a<n + l, n>l. 
The Jumarie fractional derivative has the following property: 

• The ath derivative of a constant is zero. 

• Assume that < a < 1, then the Laplace transform 

of /(") is 

£{/(«) (t)} = s a £{f(t)}-s a - l f(0). 

• (g(t)f(t)) {a) =g {a Hm)+g(t)f( a \t), 0<a<l. 

Example 1. Let /(f) = f?. Then f a \x) = r(y+ l)r~ l (y+ 1 - a)f^ a , where < 
a < 1 and y > 0. 

Example 2. The solution of the fractional differential equation 

x^ a \t)=c, x(Q) = xq, c — constant, 

is 

with the notation a! := T(l + a). 

The integral with respect to (dt) a is defined as the solution of the fractional 
differential equation 

dy = f(x){dx) a , x>0, y(0)=y , 0<a<l (1) 

which is provided by the following result: 

Lemma 1. Let f(t) denote a continuous function. The solution of the equation (Q3 
is defined by the equality 

f f{x){dx) a = a(\t- x) a - l f{x)dx, < a < 1. 
Jo Jo 

Example 3. Let /(f) = 1. Then f {dx) a = t a , < a < 1. 
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Example 4. The solution of the fractional differential equation 

=/(*), x(0)=x 

is 

x(t)= X0 + r- l (a) [\t-T) a - l f(T)dT. 
Jo 

One can easily generalize the previous definitions and results for functions with a 
domain [a,b]: 

and 

f f{T)(dTf = aj\t-T) a - l f{T)dT. 
J a J a 

For the discussion to follow, we will need the following formula of integration by 
parts: 

/ ' u {a \t)v(t)(dt) a = a\[u{t)v{t)] b a - [ u{t)v {a \t){dt) a , (2) 

Ja Ja 

where a! := T(l + a). 



3 Main Results 

Let us consider the functional defined by 

Jf(y) = [ b L(x,y(x),/ a \x),y(<*),y(b))(dx) a , 

J a 

where L(-, •) £ C 1 ([«,£>] x K 4 ;R) andx — > d^Llt) has continuous a -derivative. 
The fractional problem of the calculus of variations under consideration has the 
form 

J{y) — >extr 

(y(a)=y a ), (y(b)=y b ) (3) 
y(-)ec°. 

Using parentheses around the end-point conditions means that the conditions may 
or may not be present. 

Along the work we denote by <9,L, ; = 1,...,5, the partial derivative of function 

L(v, v, •) with respect to its ;th argument. 

The following lemma will be needed in the next subsection. 

Lemma 2. Let g be a continuous function and assume that 
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g{x)h{x) (dx) a = 



for every continuous function h satisfying h(a) — h(b) = 0. Then g = 0. 

Proof. Can be done in a similar way as the proof of the standard fundamental lemma 
of the calculus of variations (see, e.g., l40l ). 



3.1 Necessary Conditions 

Next theorem gives necessary optimality conditions for the problem (0). 

Theorem 1. Let y be an extremizer to problem (0. Then, y satisfies the fractional 
Euler-Lagrange equation 

d 2 L{x,y{x)J a \x),y(a),y{b)) = ^d 3 L(x,y(x),y^ (x),y(a),y(b)) (4) 
for all x G [a,b]. Moreover, ify(a) is not specified, then 

^d 4 L(xMx)J a \x)Ma)Mb))(dx) a ^ald 3 L(aMa)J a Ha),y(a),y(b)) (5) 
ify(b) is not specified, then 

[ b d 5 L(x,y(x),y {a \x),y(a),y(b))(dx) a = -a\dsL{b,y{b),y( a \b),y{a),y{b)). 

J a 

(6) 

Proof. Suppose that y is an extremizer of ^ and consider the value of J? at a 
nearby function y — y + eh, where e 6 R is a small parameter and h is an arbitrary 
continuous function. We do not require h(a) = or h(b) = in case y(a) or y(b), 
respectively, is free (it is possible that both are free). Let /(e) = ^ (y + eh). Then a 
necessary condition for y to be an extremizer is given by /(0) = 0. Hence, 



b r 



d 2 L(-)h(x) + d 3 L(-)h < - a) (x) + d 4 L(-)h(a)+d 5 L(-)h(b) (dx) a =0, (7) 



where (•) = ^x,y(jc),y^ a ' (x),y(a),y(b) \ . Using integration by parts (O to the second 
term we get 



b r d 2 L(-) - ^ W)] ( dx ) a + a\d,L{-)\ x=b h(b) - ald 3 L(-)\ x=a h(a) 



d 4 L(-)h(a) + d 5 L(-)h(b) (dx) a =0. (8) 
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We first consider functions h such that h(a) = h(b) = 0. Then, by the Lemma[2]we 
deduce that 

W) = £ W) 

for all x G [a,b]. Therefore, in order for y to be an extremizer to the problem ([3), 
y must be a solution of the fractional Euler-Lagrange equation ©. But if y is a 
solution of |@), the first integral in expression dH) vanishes, and then condition © 
takes the form 

h{b) [ [ b d 5 L(-)(dx) a + ald 3 L{-)\ x=b \ +h(a) [ [ b d 4 L(-)dx - a!<? 3 £(-)|.v=J = 0. 

If y(a) = y a and = are given in the formulation of problem ([3]), then the latter 
equation is trivially satisfied since h(a) = h(b) = 0. When y(b) is free, then equa- 
tion (O holds, when y(a) is free, then © holds, since h(a) or h{b) is, respectively, 
arbitrary. 

In the case L does not depend on y(a) and y(b), by Theorem [T] we obtain the 
following result. 

Corollary 1. ^ Theorem 1 ] Let y be an extremizer to problem 
Jf(y)= [ b L(x,y(x),/ a \x))(dxr 

J a 



extr. 



Then, y satisfies the fractional Euler-Lagrange equation 

d 2 L{x iy {x)^ a \x)) = ^- di L(x,y(x),y {a) (x)) 
for all x £ [a,b]. Moreover, ify(a) is not specified, then 

d 3 L(a,y(a),/ a \a))=0 

ify(b) is not specified, then 

d 3 L(b,y(b)J a \b))=0 

Observe that if a goes to 1, then the operators 4^, (dx) a could be replaced with 
jg: and dx. Thus, in this case we obtain the corresponding result in the classical 
context of the calculus of variations (see ||29l Corollary 1], |[T2l Theorem 2.1]). 

Corollary 2. 4291 Corollary 1 ] Ify is a local extremizer for 

rb 

J{y) = \ L (x,y(x),y (x),y(a),y(b))dx — > extr 
(y(a)=y a ), (y(b)=y b ), 

then 
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dx 



d^L(x,y{x),y'{x),y{a),y{b)) = d 2 L(x,y(x),y' '(x),y{a),y(b)) 



for all x 6 [a,b]. Moreover, ify{a) is free, then 

d3L(a,y(a),y'(a),y(a),y{b)) = / d 5 L(x,y(x),y (x),y(a),y(b))dx; 

J a 

and if y(b) is free, then 

fb 

d 3 L(b,y(b),y'(b),y(a),y(b)) = - / d 6 L(x,y(x),y' (x),y(a),y(b))dx. 

J a 



3.2 Sufficient Conditions 

In this section we prove sufficient conditions for optimality. Similarly to what hap- 
pens in the classical calculus of variations, some conditions of convexity (concavity) 
are in order. 

Definition 2. Given a function L, we say that L(x,y,z,t,u) is jointly convex (con- 
cave) in (y,z,t,u), if <9,L , i = 2, . . . ,5, exist and are continuous and verify the fol- 
lowing condition: 

L(x,y+yi,z+Zi,t+h,u + u{) —L(x,y,z,t,u) 

> (<)d 2 L(-)y 1 + d 3 L(-)zi+d 4 L(-)t 1 +d 5 L(-)u 1 , 

where (•) = (x,y,z,t,u), for all (x,y,z,t,u), {x,y + y\,z + zi,t + h,u + u\) e [a,b] x 



Theorem 2. Let L(x,y,z,t,u) be a jointly convex (concave) in (y,z,t,u). Ifyo satis- 
fies conditions I®-©, then yo is a global minimize r ( maximizer) to problem (O. 

Proof. We shall give the proof for the convex case. Since L is jointly convex in 
(y,z,t,u) for any continuous function yo + h, we have 

fiyo+h) - /(yo) = I \L(xMx)+Kx)Jo\x)+h (a) {x), yo {a)+h{a), 
y (b) +h(b))-L(x,y (x),y { a) (x),y (a),y (b))\ (dx) a 



> 



d 2 L(-)h(x) + d 3 L(-)h( a \x) + d 4 L(-)h(a)+d 5 L(-)h(b)]{dx) a 



where (•) = (x,yo(x) ,y^ a \x) ,yo(a) ,yo{b)j . We can now proceed analogously to the 
proof of TheoremQ] As the result we get 
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J(yo + h) - J(yo) > [ [d 2 L(-) - J^ 3 L(-)] (dx) a 
h(b) \ [ d 5 L(-)(dx) a + a\d 3 L(-)\ x=b \ +h(a) [ t d*L(-)dx - a\d 3 L(- 



= o, 



since yo satisfies conditions (IUi-©. Therefore, we obtain ^ (yo + h) > ^ (yo 



3.3 Examples 

We shall provide examples in order to illustrate our main results. 
Example 5. Consider the following problem 



r(a + l 



-2jc°y 



+(y(0)-l) 2 + (y(l)-2) 2 } (dx) a 
The Euler-Lagrange equation associated to this problem is 



extr. 



dx° 



r(a + l) 



(v 



a)\2 



2xV 



2x° 



r(a + l) 



2x a 



0. (9) 



Let y = x a + b, where Since = r(a + 1), it follows that y is a solution 

of ©. In order to determine b we use the generalized natural boundary conditions 
(O-©, which can be written for this problem as 

[\y(0)-l)(dx) a = 0, 
Jo 

[\y(l)-2)(dx) a = 0. 
Jo 

Hence, y = x a + 1 is a candidate solution. We remark that the y is not differentiable 
in [0,1]. 



Example 6. Consider the following problem 

f 



(y (a) W) 2 + 7/(0)+A(y(l)-l) 2 ] (dx) a 



(10) 



where y, A € K + . For this problem, the fractional Euler-Lagrange equation and the 
generalized natural boundary conditions (see TheoremQ} are given, respectively, as 



2|>M-0, 



(11) 



Fractional variational calculus for non-differentiable functions 9 

-1 



/ Yy(0)(dx) a = aly^(0), 
Jo 

f\(y(l) -l)(dx) a = -a\yW(l). 
Jo 



(12) 
(13) 



Solving equations (fTTb — (TT3T> we obtain that 



yXal (a\) 2 X 



7 A + (a!) 2 (A + 7)' rA + (a!) 2 (A + r) 

is a candidate for minimizer. Observe that problem ( TTOb satisfies assumptions of 
Theorem|2] Therefore y is a global minimizer to this problem. We note that when a 
goes to 1 problem ( TTOb tends to 



J^(y) = / o ' [(/M) 2 + w 2 (o) +A( y (i)-i) 2 ]^- 



mm. 



with the solution 



yX+X + y yX + X + y 



4 Conclusion 



In recent years fractional calculus has played an important role in various fields 
such as mechanics, electricity, chemistry, biology, economics, modeling, identifi- 
cation, control theory and signal processing (see, e.g., ll3l l9l [131 l25l l26l [33l l37l ). 
The fractional operators are non-local, therefore they are suitable for constructing 
models possessing memory. This gives several possible applications of the FCV, 
e.g., in describing non-local properties of physical systems in mechanics (see, e.g., 
|[T0ll24l[35ll ) or electrodynamics (see, e.g., |[8l [39l ). The Jumarie fractional deriva- 
tive is quite suitable to describe dynamics evolving in a space which exhibit coarse- 
grained phenomenon. When the point in this space is not infinitely thin but rather 
a thickness, then it would be better to replace dx by (dx) a , < a < 1, where a 
characterizes the grade of the phenomenon. The fractal feature of the space is trans- 
ported on time, and so both space and time are fractal. Thus, the increment of time 
of the dynamics of the system is not dx but (dx) a . In this note we generalize some 
previous results of the FCV (which are defined in terms of the Jumarie fractional 
derivatives and integrals) by proving optimality conditions for problems of FCV 
with a Lagrangian density depending on the free end-points. The advantage of us- 
ing the Jumarie fractional derivative lies in the fact that this derivative is defined for 
continuous functions, non-differentiable (see, Example |5J. Note that the integrand 
in problem (0 depends upon the a priori unknown final values y(a) and y(b). The 
present paper indicates how such problems may be solved. 
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